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Abstract 

Q\ The light quark-antiquark scattering Green's function is considered near a meson resonance peak. 

Q\ The Bethe-Salpeter equation is used to write formal expressions for the resonance width/mass ratio. 

t-h Arguments are made concerning to what extent this ratio can be calculated perturbatively, and an 

p | upper bound is placed on the growth of this ratio as a function of radial excitation. Certain mesons 

1) and their radial excitations are considered, as well as the more general issue of classifying mesons in 

C/} the quark model. 

^ 1. Introduction 

OS One interesting feature of QCD is the existence of a large number of meson and baryon resonances. 

Unfortunately confinement, an essential characteristic of bound states in QCD, is not completely 
^ understood and most quantitative information about these bound state resonances is deduced from 
symmetry arguments such as SU(3) flavor symmetry, chiral symmetry, and heavy quark 
f^h symmetry, or from empirical models such as quark-antiquark potential models or flux-tube 
*p models. 
Oh 

The following analysis is an attempt to deduce some consequences of perturbative QCD for light 
quark-antiquark meson resonances. The Bethe-Salpeter equation, written to lowest order in the 
perturbative expansion, is used to extract information about quark-antiquark scattering near 
resonance peaks. Perturbative approximations to the full Bethe-Salpeter equation miss important 
confinement effects, and a straightforward perturbative approximation is avoided. Instead, 
perturbation theory is employed on the expectation values of operators interposed between exact 
wavef unctions. The assumption is made, and motivated later, that the expectation values of non- 
perturbative operators associated with confinement can be neglected. 

2. Quark-Antiquark Scattering and the Bethe-Salpeter Equation 

The Green's function of interest is 

Gr(/W,«M1/>-//) = (o^ (2.1) 
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where the Greek superscripts are Dirac indices, the contracted Latin indices show the color 
singlets, and the subscripts 1, 2 refer to quark flavor. Only connected diagrams contributing to 
this process will be considered, as they should be sufficient to describe properties of the meson 
resonance peak. 

As this Green's function is gauge dependent, many of the quantities derived from it will also be 
gauge dependent, including the Bethe-Salpeter wavefunctions. However, physical quantities such 
as masses and widths will, as usual, be gauge independent. Unfortunately, the rather useful 
process of checking the gauge invariance of physical quantities explicitly will not be possible, 
because wavefunctions are not calculated. 



Two particle irreducible (2PI) diagrams are defined as those that do not contain a purely quark- 
antiquark intermediate state. By convention, the external quark and antiquark lines of 2PI 
diagrams are removed. Below is a diagrammatic recurrence relation, known as the Bethe-Salpeter 
equation, for the full connected Green's function written in terms of 2PI diagrams and full quark 
and antiquark propagators [1]. 
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By using a matrix shorthand in which matrix multiplication signifies contracting Dirac indices and 
integrating over momenta, the recurrence relation can be written succinctly and solved. 

G = K + KIG 

g={k-'-ij' 



(2.2) 
(2.3) 



G = full connected Green' s function 

K = full quark and antiquark propagators 

/ = two particle irreducible Green' s function 

Let 0. be the eigenfunctions of G, and A, be the respective eigenvalues. 

G0. = A.0. 

It should be noted that the matrix G will, in general, not be hermitian. Let 0., represented as a row 
vector, be the eigenfunctions of G T , and A be the respective eigenvalues. 



(2.4) 



0. G = A. 



(2.5) 



As any matrix and its transpose share the same characteristic equation, they share identical 
eigenvalues. Let M be the subspace spanned by the nondegenerate eigenfunctions of G. In this 
subspace, the indexing of both sets of eigenfunctions can thus be aligned by matching eigenvalues. 
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index/: Gfr = Xrf,, <j) t G = (2.6) 

Since X,fa = j>G<j> t = X k fa fa = cd ik . (2.7) 

Let 5(5) be the subspace of M spanned by <j) t ( 0.) for c nonzero. Using orthogonality relations, the 
eigenfunctions can be used to write a bilinear form, G' , equivalent to G on these restricted 
subspaces. 

G' = £A,(#,f#, (2.8) 

fGg = fG'g VfsS, gsS (2.9) 
By the Bethe-Salpeter relation derived above, 

A; 1 = (fa^ = (00) ^.(/r 1 - 7)0 = (^.j'^jr 1 ^ - 0/0) (2. 10) 

ff 3 !^-^ (2 - n) 

3. The Resonance Peak 

By momentum conservation, the eigenfunctions described above can be chosen to be states of 
definite total momentum. To avoid confusion, this total momentum will be written as a subscript 
to make explicit its role as an index and to distinguish it from the relative momentum in whose 
Hilbert space the eigenfunctions reside. 



G'(p;-cf,q) = Y^ 



(3.1) 



where f,(p 2 ) = j] ctq ctq\\ p {q')K \p; - q', q^ J-q) - $Jq)l{p; - q, q)<P, J-q)) (3.2) 

Assuming that f.(p 2 ) is well-defined and that \f,(p 2 )\ has a unique minimum, the location of that 
minimum will be denoted 

fa \flfa)\ = min{\f{p 2 )\} (3.3) 
If |./](# 2 min )| is small then G'(p; - q, q) near p = p 2 min is dominated by 

[fAp^YhJ-rtkb) ( 3 - 4 ) 

By adjusting the eigenfunction's phase, the behavior of f t (p 2 ) near p 2 min can be written as 

f{p 2 )~a(p 2 -p 2 J + ibp 2 Amn a,b real (3.5) 

The interpretation of these results is that the square of the mass of the resonance is p 2 min . The 
width/mass ratio of the resonance is b I a. 



4. Perturbative Calculations 

This section applies the formalism of the preceding sections to meson resonances. Several 
assumptions and approximations are made. These assumptions and approximations are motivated 
by qualitative physical arguments, but the actual size of the errors incurred are not known. 

Even for highly excited meson resonances, where energies are sufficient to suggest that the quark- 
antiquark interaction is predominantly in the asymptotically free regime, confinement must still 
modify the behavior of the wavefunction in the region where the quark and antiquark are widely 
separated. The present analysis adopts the view that for such highly excited meson resonances, 
confinement can be effectively modeled as a large distance boundary condition on the 
wavefunction. 

If this confinement boundary condition is sufficiently abrupt and severe, a high wall for instance, 
the wavefunction will not penetrate substantially into regions in which confinement is important. A 
simple analogy can be made using non-relativistic quantum mechanics in one dimension. If one 
starts with a smooth potential V and adds to it a pair of infinite walls at locations a and -a, the 
eigenstates of the new potential, V, will vanish at a and -a. For the purposes of calculating the 
expectation value of V, the walls can be neglected, as, in some sense, information about the walls 
are contained in the wavefunctions themselves. It will be assumed that something similar happens 
in the case of highly excited mesons, and the expectation values of non-perturbative, confinement- 
producing operators are neglected. 

All perturbative calculations are done in Landau gauge and dimensionally regulated in the MS 
prescription. Light quark masses are neglected although, to some extent, they make a reprise in 
phase space corrections described later. The running strong coupling constant, renormalized at the 
mass of the resonance under consideration, m, is used to minimize higher loop corrections. In 
general, diagrams are kept to lowest possible order. Free quark and antiquark propagators are 
used in the matrix K. A freely propagating t-channel gluon is the lowest order contribution to the 
dispersive part of the matrix /. The absorptive part of the matrix /, however, is a more subtle 
matter and gets its lowest order contribution from a higher order diagram. 

The absorptive part of any process is associated with intermediate states situated on their physical 
mass shells. As the physical states of the underlying free theory are quite different from the 
physical states of fully-interacting QCD, naive lowest order perturbative calculations of the 
absorptive part of a process can be misleading. For instance, in calculating the propagator for 
meson resonances, the lowest lying possible physical intermediate states in the fully-interacting 
theory are two meson states. In order for the intermediate two meson state to be produced, an 
extra quark- antiquark pair must first be produced since a two meson state has two valence quark- 
antiquark pairs. The simplest approach then, and the one followed here, is that the lowest order 
contribution to the absorptive part of a strongly interacting process is given by the lowest order 
diagram(s) with the correct valence particle structure. This will be referred to as valence 
correspondence. In general, three body or higher body decays occur less frequently, and only 
decays into two mesons are treated in this analysis. . Below are the lowest order diagrams and 
calculations for the expectation values of equation (3.2). 
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(Full) pull) 
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Lowest Order 
(absorptive) 



A 



o 



Evaluating these diagrams, 

K(p 2 ) - ^>,K-% p « i ■ j d 4 q (-tr(<t>. p (q)[p - q]^ p (-q)[p + q])) (phase chosen to be real) 
l{p 2 ) = fob_ f « 4 j| A A j Re[ g 2 (m 2 ) • M p (q, q')} + i • • Imflnf ~^~^7~ /£ ll • M p (q, q') 



where M p (q, q') = i 



($-$0+/e ((4_4') 2 + ; e } 
At this point, it is convenient to define some auxiliary functions, 

Xl (p 2 ) = Rejo lo } = Z +(p 2 ) + ^(z? 2 ) 

W 2 ) 



fr[^)^.-,(-«)/v] 



^).lm{^} = i .„ / .^)- 



g 2 (m 2 )-0 



where ^(p 2 ) = 4 • RejjJ ctq ctq 
Equation (3.2) can now be written as 

f{ P 2 ) = k{ p 2 ) - [ %l (p 2 ) = x ;(p 2 )+ x;(p 2 )] + 



rrt 



'\ 2 \ 



■M(q,q') 



„, 2, • S 2 {m 2 ) 2 
x(p) = i-n f -^^-X i (p) 



(II) 



(HI) 



(4.1) 
(4.2) 

(4.3) 



(4.4) 
(4.5) 

(4.6) 



(4.7) 
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n f refers to the number of possible flavors of light quark loops in the absorptive diagram shown in 
figure (III). The simplest approach is to set n f =3. It may be useful, however, to be more careful. 
Although the masses of light quarks have been neglected in the above calculations, the mass of the 
light quark does become important near a kinematic threshold. One example is the p(770) , which 
can decay into two pions but is kinematically forbidden to decay into two kaons. To satisfy 
valence correspondence (VC), the s-quark loop should not be included in counting n f for the 
p(770), since there are no kinematically allowed two meson states with the correct valence 
properties. In order to handle such kinematic effects systematically, the following formula, 
motivated by the phase space factors for two body decays, will be used to evaluate n f . 



n, ~ J max 

n _ [wo meson states 

flavor ofqq (nmc , ni| mj) obi:yillg vc 



[m - (m, + m 2 J )(m 2 - (m, - m 2 ) 2 j 



m 



(4.8) 



It is useful to further analyze K^p 2 ). By charge conjugation and SU(3) flavor symmetry, terms 
linear in q in equation (4.1) do not contribute. 

K{p 2 )~^{p 2 )+^{p 2 ) (4.9) 

where <(p 2 ) = i- J ctq tr[-<j>. p {q) p ftj-q) p ] (4.10) 

and Kf^^i-j^trf-i^^^H)?] (4.11) 

The total spin of quark-antiquark mesons arise from the intrinsic spin of the quark and antiquark as 
well as their orbital angular momentum. It is known that the intrinsic spin of the quark-antiquark 
pair is either spin one or zero. This fact is used to simplify equation (4.10). 

<t>,, P = K P ■ /. + <pZ • /.. + <pU ■ i + r • r, • r 5 (4.12) 

o =<T. ./.+<? .. /.-<p\-l + r -7 -7 (4.13) 

K:{p 2 ) = p 2 -d,{p 2 ) (4.14) 

where d{p 2 ) = \ d 4 q [^{q)^{-q)+ ¥MW,{-q)+ f(?KH)+ ^" (?)<(-?)] (4-15) 

The next step is to write down an expression for f,(p 2 ) near p 2 = m 2 Using the orthogonality 
condition derived in equation (2.7) and assuming completeness of the Bethe-Salpeter 
eigenfunctions at p 2 = m 2 , f(p 2 ) can be expanded in terms of eigenfunctions at p 2 = m 2 . The result, 
to first order in p 2 -m 2 , is simply 

f : (p 2 ) - JJ ctq d^iq^K'ip; - q\ q)fr{-q) - j>{q)l{p; - q, q)^{-q)) (4.16) 

where ^ = ^^ and ^, = ^__ m! (4.17) 

This can be rewritten in terms of the various auxiliary functions defined above, 

/, {p 2 ) - K ■ ^(m 2 ) + K?{m 2 )-%:{m 2 ) - ^(m, 2 ) + n(m 2 ) (4.18) 

m. 

<(m i 2 ) + K*(m 2 ) - %;(m 2 ) - X-{m 2 ) - (4.19) 
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The resonance mass/width ratio, r , is given by 



r = 



<(m 2 ') 



Using equations (4.8), (4.14), and (4.19), 
g2 ( m -} rim 1 ) 



r ~ 



max 



_ Iwo meson states 
flavor ot qq (mas ses m,, ra,) obeying VC 



{ni - (m, + m 2 J )(m 2 - (m, - m 2 )" j 



m 



m = 



K p . (m. 2 ) _ - Kf (m. 2 ) + X-{m : 2 ) + %: (m. 2 ) 



(m. 2 ) 



<f (m. 2 ) 



(4.20) 



(4.21) 
(4.22) 



5. Radial Excitations 



It is convenient to define the normalized quantities, 



K-f(m, 2 ) 
d i (m,. 2 ) 



The remainder of this analysis will be concerned with the behavior of these three quantities as a 
function of radial excitation number, n, with all other parameters fixed. It is assumed that the mass 
of any sequence of radial excitations of resonances is, for sufficiently large n, well approximated 
by a smooth, monotonically increasing function of n. Since -i^{m 1 ) + x*{m 1 ) + xj{rn) = rn, it is 
also assumed that each of -Kf(m 2 ), X-{ m ~)> and fcirf) can ^ e approximated by smooth functions 
of n. Furthermore, it is naively assumed that for large n these functions are well-behaved and do 
not oscillate rapidly nor with wide amplitudes (e.g., power law or logarithmic). 

The following definition will be useful, 

k. = 



^4tT 



7. max 

n _ two meson states 

flavor of qq (mlsses ^ ^ obeying yc 



{rri — (m, + m 2 ) )(m 2 - [m { - m 2 ) j 



m 



-K*{m) + x+{m) + x;{ni) 



(5.1) 
(5.2) 



(5.3) 



(5.4) 



As n becomes large, one possibility is that fairf) dominates -Kf(m. 2 ) and x\ m T)- Sequences of 
radial excitations of meson resonances with this property, if there are any, will be denoted as type 
1. Type 1 sequences have the property that for sufficiently large n, & — > 1 . 

A second possibility is that X*(rn) is dominated by either -Kf(/rz 2 ) or X-irf)- These resonance 
sequences will be denoted type 2 and have the property that for sufficiently large n, & — > 0. 
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The third possibility is that %*{ m ^) shares dominance with one or more of the other two terms, i.e., 
has the same large n behavior. Resonance sequences with this property will be denoted type 3. 
For large n, type 3 sequences will have & — > & , not necessarily or 1. 

For all three types of sequences, k. is bounded above. This is simply the statement that the 
denominator in equation (5.4) is not fine-tuned to be arbitrarily smaller than the numerator. The 
bound on k. gives a bound on the growth of r. . 



3k >0 s.t. r<k- 



"24tT 



7. max 

n _ two meson states 

flavor of qq (masses ^ ^ obeying yc 



[m - (m t + m^[m — (n\ - mj j 



m 



(5-5) 



6. Experimental Data 

There are a number of mesons for which a reasonable quark model interpretation exists, and the 
initial entries of some radial excitation sequences can be examined. The data for these is presented 
below [21. In the calculation of L, A— is taken to be 250 MeV. 

* Mi 





7(G)/ TP(C)\ 
1 \ J ) 


25+1 j of V m nHp 1 ^ 


r i 




p(770) 


1 + (1~) 


\\ 


.197 + .0016 


1.02 + .01 


p(1450) 


1 + (1~) 


2\ 


.212+. 041 


1.191.23 


JT(892) 


Hi") 


\% 


.0558+. 0009 


.471+. 008 


iT(1410) 


Hi") 


2\ 


.161+. 016 


.996+. 099 


0(1020) 


0(1-) 


\\ 


.00434+. 00006 


.103+. 001 


0(1680) 


0(1-) 


2\ 


.0893+. 0298 


.616+. 206 


©(782) 


0(1-) 


\% 


.0108+. 0001 


(?) no 2-body decays 


(0(1390) 


0(1-) 


2\ 


.123+. 042 


.880+. 304 


/ 2 (1270) 


+ (2 ++ ) 


VP 2 


.146+. 016 


.793+. 087 


/ 2 (1810) 


+ (2 ++ ) 


2T 2 


.214+. 010 


1.301.06 


£(1525) 


+ (2 ++ ) 


VP 2 


.0498+. 0066 


.350+. 046 


£(2010) 


+ (2 ++ ) 


2T 2 


.100+. 032 


.697 + . 32 


£(1430) 


H2 + ) 


VP 2 


.0688+. 0016 


.427+. 010 


£(1980) 


H2 + ) 


2T 2 


.197+. 020 


1.291.13 


X 


1"(0" + ) 


VS 


stable 


°o in chiral limit 


7r(1300) 


1"(0" + ) 


2\ 


.338+. 062 


3. 81+. 70 


tt(1770) 


1"(0" + ) 


3'5 


.175+. 028 


1.341.21 


K 


HO") 


1'5 


stable 


°o in chiral limit 


K(U60) 


HO") 


2\ 


-.178 


-1.96 


^(1830) 


HO") 


3S 


-.137 


~1.11 



The chiral limit divergence of k for % and K is obtained by noticing that the denominator in 
equation (5.8) vanishes as m] or nf K , whereas the numerator, in general, need not vanish. It may 



be argued, qualitatively, that the large values of k. for the 2 S pseudoscalars is a consequence of 
the relative smoothness of k. as a function of radial excitation. 

Higher excitations of mesons will have to be examined in order to check if k does actually appear 
to converge and, if so, what k converges to in order to classify the excitation sequence as type 1, 
2, or 3. 



7. Remarks 

The main theoretical result of this analysis is that the width of mesons in a radial excitation 
sequence does not grow arbitrarily large with respect to meson mass, and the upper bound on the 
growth of the width is given in equation (5.5). It is worth noting that this result should be valid, 
with slight changes of parameters, for confining QCD-like theories such as technicolor. 

If in the future a radial excitation sequence of reasonable length is identified, classifying this 
sequence as type 1, 2, or 3 may be useful in providing some information about the distribution of 
the wavefunction in momentum space. For immediate purposes, however, it may prove useful to 
employ the converse and assume that k is convergent in order to check whether a given meson is a 
good candidate for the next radial excitation of a given sequence. In this regard, one might have a 
potentially useful tool to help classify mesons in the quark model. Unlike SU(3) flavor or nonet 
symmetry which allows "horizontal" comparisons of mesons within an approximately degenerate 
multiplet, this technique would be a "vertical" comparison of radial excitations of mesons with 
sequentially increasing mass. 
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